Determination and prediction of the fretting crack initiation: introduction of the (P,Q,N) 
representation and definition of a variable process volume 
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In this work, the crack nucleation under fretting loading is investigated experimentally with a damage tolerant 
2024 aluminium alloy. A new method is introduced to determine its condition with respect to all loading 
parameters including the number of fretting cycles. Further work deals with the prediction of this threshold 
using the Smith- Watson-Topper criterion. New developments are presented, in particular a process volume of 
variable size is introduced in the computations of the fretting crack initiation. 
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I. INTRODUCTION 

Fretting damage has been recognized as a problem in sev- 
eral industrial applications for years now. The first studies 
defined the concept of fretting [1], while further work intro- 
duced more rationalization with precise concepts such as fret- 
ting cycles and fretting maps |2|,l3|]. Depending on the contact 
parameters and especially the displacement amplitude, it has 
been shown that two main regimes can be defined. First, In 
the large displacement range, the full sliding regime exists, 
which can induce wear of the surfaces in contact. Second, In 
the low displacement amplitude range, the contact causes a 
partial slip condition and severe stress gradients at the contact 
border. This can induce very rapid crack nucleation compared 
to classical fatigue testing. In this paper we focus on the lat- 
ter regime and particularly on the crack nucleation condition. 
This has been recognized as a critical issue for industry in the 
past ten years. 

Thanks to recent work, different methods are now avail- 
able to model the propagation of fretting fatigue cracks using 
for example the "Crack analogue" model, fracture mechanics 
0, Hi and also weight functions jfj. But due to the possi- 
ble drastic reduction of the fatigue limit induced by fretting 
10], a great part of the work has been devoted to determin- 
ing the fretting crack initiation conditions in various types of 
materials |[8il9[ [Toll . This process is quite complex due to high 
stress and strain gradients, oxidation phenomena and other tri- 
bological phenomenon (TTS formation [11]), surface defor- 
mation and wear, and debris formation. Until now the main 
method used to describe quantitatively crack nucleation deals 
only with the mechanical state under the fretting conditions 
and relies on computing a multiaxial fatigue criterion (this 
was first attempted by Petiot et al. Il2tl ). This can be done 
either by using analytical solutions of the stress/strain state 
under fretting contact (as it can be found in 11311 ). or more re- 
cently by numerous author with finite elements iTToL 0, fl5ll . 
The latter method is more constrained but allows one to test 
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more general geometries and more complex loading condi- 
tions. Due to the very strong gradients located at the contact 
border, one must use a process volume approach, which con- 
sists of averaging the stress/strain state in a micro-volume of 
material prior applying the multiaxial criterion. This method 
has popularized and has been investigated with various fatigue 
criteria |@, EH El • These investigations have been capable of 
predicting, in some cases, the crack nucleation threshold ; but, 
the significance of the process volume is not yet fully under- 
stood, and as a consequence cannot be truly estimated with- 
out doing experiments. For completeness, it must be pointed 
out that some alternative methods exists such as the asymp- 
totic stress intensity analysis 01711 . In this kind of approach, 
the stress singularity induced by the contact configuration is 
captured by asymptotic analysis in order to deduce a stress 
intensity factor at the contact edge. 

In this paper we restrict the analysis to the prediction of the 
crack nucleation condition in the fretting wear configuration 
by the Smith- Watson-Topper criterion. We present an new vi- 
sion of the process volume and apply the method to a complete 
experiment series with various fretting wear conditions. This 
work attempts to clarify the physical meaning of the process 
volume and to set up a systematical way of characterizing the 
crack nucleation condition with respect to parameters relevant 
for the industrial applications (contact pressure, shear stress 
amplitude and number of cycles). 



II. EXPERIMENTAL WORK 

Part of the experimental work has been presented else- 
where. The reader is referred to [16] for more details on 
the experimental setup, though general information is sum- 
marized here. A general fretting wear apparatus is used in 
cylinder/plane configuration and under partial slip regime. A 
normal force P is applied on the counter-body to maintain the 
surface in contact causing an elliptic pressure p(x) over the 
contact zone |i| < a ; a cyclic relative displacement is im- 
posed to the interface leading to a classical shear stress q(x) 
which exhibits its maximum at the stick zone limit x = ±c (cf. 
fig-ID' The flat specimen being maintained in a fixed position, 
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this displacement gives rise to a cyclic shear force 1 of magni- 
tude Q which is measured during the test via a force sensor. 
Regarding the materials, an aluminium/aluminium contact is 
studied, with plane samples in a 2024T35 1 alloy (see table Q] 
for mechanical properties) and a cylindrical counter body of 
49 mm radius made of a A17075T6 alloy. 

TABLE I: Mechanical properties of the studied A12024 



E (GPa) 



0o.2% (MPa) a d (MPa) 



72 0.33 



325 



140 




FIG. 1: Schematic of the fretting wear cylinder/plane configuration. 

The fretting crack nucleation has been investigated in the 
following way (cf. fig. 01: 

• fretting test is conducted on the sample for chosen load- 
ing conditions (fig.|2^). 



newly created face (fig. |2j>,c). 

chemical etching can be performed after a first opti- 
cal micrograph observation in order to avoid a possible 
blurring of the crack by the polishing process (fig. 

depending on the previous investigations, new fretting 
conditions are determined to refine or confirm the re- 
sults. 




FIG. 2: Experimental method to investigated cracking after a fretting 
test ; see the text for details. 



A previous study (with the same geometry and material) 
conducted to 50000 fretting cycles was described in llrjll . 
According to these experiments, the crack nucleation condi- 
tion was found to be mainly independent of the normal load 
(see fig. [8] where the experimental results have been replot- 
ted), leading to the crack nucleation threshold in terms of the 
tangential load Q t h — 240 N/mm. 

Though this value defines a boundary between safe and 
a crack risk zones, the number of cycle needed to be more 
deeply investigated. Indeed a rapid calculation of the stress 
state at the contact border does show that the contact is sub- 
jected to loading conditions equivalent to low cycle fatigue. 
We consider the simple analysis of the two-dimensional plane 
strain cylinder/flat contact with the Hertz theory from which 
we can write the peak contact pressure pq: 



Po = 



2P ( PE* 



1/2 



(1) 



R being the radius of the cylinder and E* the effective 
Young modulus. 

The stress field is known in partial slip conditions from the 
analytical solutions (see for instance [13]). At the contact bor- 
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matches the principal stress state (a±, 02, 03) which is biaxial 
(the contact pressure vanishes): 



(2) 



with fjb the friction coefficient in partial slip condition (fi 
was determined to be close to 1,1 in our case Biol ). The equiv- 
alent Von Mises stress a e can be easily expressed: 



1 
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combining with equation [2] 



a e (x = a,y = 0) = 2p 



2po 




fiE*Q 2 



(4) 



1 note that in fretting wear experiments, there is no bulk load imposed 



Thus the shear force leading to the conventional yield at 
0.2% of deformation cry, in the slip zone, can be written as: 
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With cry = Co. 2%' this expression gives Qy = 117 N/mm. 
It is clear that an aluminium/aluminium contact with a fric- 
tion coefficient greater than unity (like in our case) will in- 
duce very high contact stresses and it can be seen that all the 
tests carried out in this study operates under low cycle fatigue 
range (Q > Qy)- It is therefore expected that the number of 
fretting cycles must have an influence on the crack nucleation 
condition. 

To fully determine the crack nucleation condition, the ex- 
perimental study was completed as follows: for a chosen nor- 
mal pressure, additional fretting tests were conducted with 
different numbers of cycles. The initiation threshold is then 
determined for each new number of fretting cycles. The crack 
nucleation boundary is then expressed as a function of the 
three loading parameters P, Q and N. The whole experimen- 
tal procedure is described in fig. [3] 




FIG. 3: Methodology used to determine the crack nucleation condi- 
tion with respect to the ,Q,N parameters, a) tests at constant N, b) 
tests at constant P, c) whole (P,Q,N) representation. 

Here, new tests have been performed with a constant nor- 
mal load P = 400 N/mm. Two different numbers of fretting 
cycles have been investigated: 5.10 5 and 10 6 cycles. The cor- 
responding crack nucleation threshold are determined accord- 
ing to fig. [3^ ; values have been reported in tableland plotted 
on fig. ED 

TABLE II: Influence of the number of fretting cycles on the crack 
nucleation threshold (P = 400 N/mm). 



Number of cycles 


Threshold (N/mm) 


5.10 4 


240 


5.10 5 


190 


1.10 6 


170 



As expected, the number of applied fretting cycles has 
an influence on Qth- The more fretting cycles applied, the 
less the threshold value. Moreover, the evolution of Qth(N) 
shows first a rapid reduction with increasing N, and then tends 
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FIG. 4: Evolution of the critical tangential force Q t h to initiate a 
fretting crack (determined as shown in fig. [3^), with the number of 
fretting cycles (P = 400 N/mm). 

to a saturation value around 160-170 N/mm , although more 
experiment are needed to establish the exact value. In partic- 
ular, experiments with N > 10 6 cycles, which are very time 
consuming, are required. In addition, it has been shown else- 
where that with this material in this configuration, 10 6 cycle 
does correspond to a stabilized condition 111 (811 . 

Assuming a homotetic behaviour with respect to P, the 
whole crack nucleation condition can be plotted in a 3D di- 
agram f(P, Q, N) as shown in figure[5] 




2DOO ' O 



FIG. 5: Representation of the crack nucleation condition in the 
(P,Q,N) space. 

This defines completely both the safe and cracking domains 
by the boundary f(P,Q,N). The advantage is the strong 
physical meaning of these parameters, which can be extended 
to other test configurations and more generally to any struc- 
ture. In fretting fatigue tests for instance, the definition of P 
is unchanged but there is an additional load (bulk stress) in 
the material; the parameter Q must take that bulk stress into 
account and N would be the number of fatigue cycles. In a 
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more complicated case like an industrial component, first step 
will be to locate the contact area and then calculate the con- 
tact pressure P, for example by finite elements (integrating 
the normal stresses over the contact surface). Then Q can be 
defined as the maximum shear force transmitted through the 
contact during a cycle and should be estimated in the same 
way, integrating this time the shear stresses over the contact 
surface (this will, however, require a detailed contact analysis 
which must capture the partial slip contact conditions). N can 
then be identified as the number of cycles associated to the 
load responsible of the relative displacement of the surfaces. 
To conclude, this new representation aims to narrow the gap 
between ideal test configuration and actual structures encoun- 
tered in real industrial problems. 



III. PREDICTION OF THE CRACK NUCLEATION 
CONDITION 

The purpose of this section is to use a multiaxial fa- 
tigue criterion to predict the crack nucleation condition. The 
stress/strain state is computed from the analytical solutions of 
the cylinder plane contact ; these fields are averaged 2 on a cir- 
cular shaped process volume to soften the contact gradients 
and then the cracking risk is evaluated. Among the numerous 
studies on the subject, it has been shown that no criterion is 
able to predict the location of nucleation, the crack angle and 
load leading to crack initiation together. 

In this study we consider the SWT parameter despite the 
fact that it has been shown to be unable to predict the crack 
initiation angle (the investigated cracks are very often inclined 
toward the centre of the contact). The hypothesis to identify 
the crack initiation plane as the critical plane of the criterion 
is not obvious and the crack initiation process may be much 
more complicated as recently pointed out by another study 
fil^l . On the other hand the SWT criterion has a low cycle 
fatigue behaviour which has been shown to be effective for 
finite life H20TI and also was used recently be Fridrici et al. to 
tackle the same kind of problems in titanium alloys 112 ill . In 
addition, the material coefficients for this criterion are given 
by Szolwinski et al. H. 



A. Prediction of the crack nucleation boundary at 50000 
fretting cycles 

The use of the SWT criterion for fretting is now well doc- 
umented (see for example lfl5l[T6ll2lll ). Thus, only the gen- 
eral equation is reproduced here. According to this model, the 
initiation is likely to occur on the plane where the SWT pa- 
rameter is maximum. The SWT parameter T is evaluated in 
the current plane by the product between the amplitude of the 
strain e a and the maximum stress normal to this plane <r max . 



T = o max xe a = { ^-(2N) 2b ' + a' f e' f (2N) b '+ c ' (6) 

where is the fatigue strength coefficient, b' is the fatigue 
strain exponent, e'j is the fatigue ductility coefficient and c' 
is the fatigue ductility exponent. The mechanical and fatigue 
properties of the studied alloy are listed in tables |T] and Hill 
respectively. 

TABLE III: Fatigue properties of the studied A12024 (from d) 
a' f (MPa) b' e' f c' 
714 -0.078 0.166 -0.538 

This approach has been previously used to study the crack 
nucleation at 50000 cycles iflrjll . The results can be summa- 
rized as follow: 

• the crack nucleation threshold predicted by the analy- 
sis of the local analytical stresses is very far from the 
experimental result (~ 100 N/mm to be compared to 
Q th = 240 N/mm). 

• performing size effect calculations allowed for fitting 
of the average experimental result but the precise be- 
haviour (i.e. the effect of the loading pressure) cannot 
be captured. 

• the average process volume showing the best correla- 
tion was found to be correlated to the mean grain ra- 
dius (which was measured by Electron BackScattered 
Diffraction analysis giving 75 |jm). 

The physical meaning of the process volume size is still an 
open question. The most often quoted significance is the grain 
size. Indeed it is argued here that the stress/strain state induc- 
ing initiation must be sufficiently widespread to make the very 
short crack propagate in the adjacent grain. Looking at several 
studies on different materials, a large range of sizes have been 
used for the process volume. 5 pm in steel llloll . 30 pm in tita- 
nium alloy [21] and 80 |jm in aluminium alloy [16], each time 
correlated to the grain size. On the other hand, the grain size 
cannot be the only significant parameter due to other cases 
where no correlation has been found. Moreover, this param- 
eter does not take into account any mechanical quantity such 
as for example the slip amplitude which is certainly relevant 
for the crack initiation condition. 

We consider here a new approach using a variable process 
volume. This idea comes from the fact that the crack initiation 
may be strongly monitored by the stress state in the slip zone 
only. More precisely, the severity of the stress gradient located 
in the slip zone may be responsible of the crack initiation risk. 
This is illustrated by figure [6] where two distributions of the 
surface traction 3 {a xx ) are plotted for the same conditions but 



2 this process is sometimes referred to as a size effect 



3 One should note that the <j xx component is dominant in the slip zone, 
compared to the other stress values, (they vanish when x t— » a where the 
stress state becomes purely uniaxial) 
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with different normal loads (Pi < P-£). 



Pi = 240 N/mm 
a„ (MPa) Qi = 24 N/mm 




P 2 = 500 N/mm 
j„ (MPa) O2 = 240 N/mm 
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FIG. 6: Illustration of the radius of the process volume needed to 
average the stress gradient with two different normal loads Pi and 
Pi. 
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FIG. 7: Flowchart of the different steps required to predict the critical 
load Qth for each pressure level Pi ; the third step has been added to 
the conventional constant process volume approach. 



It is clear that the severity of the stress gradient can be re- 
lated to the slip zone width a — c due to the quasi-linear evolu- 
tion of a xx . Keeping the process volume size constant would 
introduce a strong effect of the pressure, which can actually 
be seen looking at the solid line prediction replotted on fig. [8] 
for the calculation with r = 80 pm. From here the radius of 
the process volume zone is no longer constant. We introduce 
a new variable 7 defined as the ratio of the process volume 
radius r and the slip zone width a — c: 



1200 



7 = 



(7) 



The value of 7 is identified once for all for the conditions 
which gave r=80|im. For the corresponding loading condi- 
tions (P = 318 N/mm, Q = 240 N/mm) one can calculate 
a = 700 pm and c = 393 |jm. This leads to the value of 
7 = 0.26. The crack nucleation threshold is then computed 
through this new approach according to the flowchart shown 
infig.0 

Prediction of the crack nucleation boundary for 50 000 fret- 
ting cycles with 7 = 0.26 is presented in fig. [8] This approach 
gives, a very good correlation with the experimental results. 
In particular, the pressure effect is well described. 

This result highlights the reliability of the approach and 
shows that the process volume size, in addition to be related 
to some microstructural characteristic of the material, must be 
linked to the mechanical state as well, such as the slip zone 
width, to be able to predict the crack nucleation precisely. 



B. Prediction of the low cycle fretting behaviour 

The same model is applied to test the experimental results 
obtained for 5.10 5 and 10 6 fretting cycles (see jffTt. The 7 pa- 
rameter is kept constant at the same value identified for 50 000 
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FIG. 8: Prediction of the crack nucleation boundary at 50 000 fretting 
cycles with the variable process volume approach. 



cycles (7 = 0.26, see jjLLLAI . The prediction is done in the 
same way by varying the number of cycles in equation|6] The 
results of these computations are gathered in figure [9] 

The experimental behaviour appears well correlated with 
the prediction of the SWT parameter. In particular one can see 
a very rapid decrease of the critical tangential force needed to 
nucleate a fretting crack, as the number of cycles increases. 
Around 10 6 cycles a plateau is reached, corresponding to en- 
durance conditions. This behaviour is consistent with another 
study on the same alloy and another 7xxx series alloy d 
where 10 6 fretting cycles are clearly identified as a stabilized 
condition. 
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FIG. 9: Prediction of the evolution of Q t h with the number of cycles, 
by application of the SWT criterion suited with a variable process 
volume size. 



IV. CONCLUSION 



Two main results have been presented in this paper. On the 
experimental point of view, a new representation of the crack 
nucleation condition is introduced through the (P, Q, N) di- 
agram, leading to a concrete and complete description of the 
material fretting resistance to initiation and this with a lim- 
ited, although still quite important, number of tests. In or- 
der to predict this initiation condition, the SWT criterion is 
used and the classical computation is extended with a variable 
process volume size. This further opens the question of the 
significance of this parameter as it appears not to be only re- 
lated to a microstructure characteristic length. The use of the 
slip zone width to determine the process volume radius clearly 
identifies a mechanical significance. Eventually the final an- 
swer may require a mix of the microstructural/mechanical be- 
haviour. This could at last be answered by extending the ap- 
proach to different materials as for instance steels and titanium 
alloys. 
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